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Abstract 


A  generalized  likelihood  ratio  test  is  known  to  be  able  to  reliably  detect  a  signal 
known  except  for  amplitude  in  incompletely  characterized  colored  non-Gaussian 
noise,  although  it  is  computationally  intensive.  A  Rao  efficient  score  test  is  proposed 
as  a  computationally  simpler  alternative.  The  Rao  test  shares  all  the  asymptotic 
properties  of  the  generalized  likelihood  ratio  test  for  large  data  records  and  small 
signal  amplitudes.  Its  detection  performance  is  asymptotically  equivalent  to  that 
obtained  for  a  similar  detector  designed  with  a  priori  knowledge  of  the  unknown 
noise  parameters.  Computer  simulations  of  the  performance  of  the  Rao  detector 
support  the  theoretical  results.  A  Rao  detector  built  with  the  knowledge  of  the 
true  form  of  the  noise  PDF  is  shown  to  significantly  outperform  a  detector  which 
assumes  the  noise  to  be  Gaussian. 


I.  Introduction 


Detection  of  a  weak  signal  in  noise  is  a  problem  of  general  interest,  having  been 
addressed  by  previous  researchers.  The  problems  studied  range  from  detection  of  a 
completely  known  signal  in  known  white  Gaussian  noise  [Van  Trees  1968]  to  detection 
of  an  unknown  signal  in  unknown  colored  non-Gaussian  noise.  An  attempt  to  solve  the 
latter  problem  has  been  made  recently  by  Kay  and  Sengupta  [1986,  3].  It  was  assumed 
that  the  signal  is  known  except  for  its  amplitude  which  can  be  positive  or  negative.  The 
problem  is  cast  as  testing  of  composite  hypotheses  and  a  generalized  likelihood  ratio 
test  (GLRT)  [Kendall  and  Stuart  1979]  is  employed  to  solve  it  a  using  parametric  rep¬ 
resentation  of  the  noise  statistics.  The  GLRT  is  found  to  be  particularly  well-suited  for 
this  problem  in  that  it  has  many  attractive  asymptotic  properties  such  as  consistency, 
unbiasedness  and  constant  false  alarm  rate  (CFAR).  It  is  also  asymptotically  optimal  in 
the  sense  that  knowledge  of  the  unknown  noise  parameters  would  not  improve  its  per¬ 
formance  provided  the  data  record  length  is  sufficiently  large.  However  the  GLRT  has 
a  serious  disadvantage,  namely,  its  computational  complexity.  It  necessiates  computing 
the  maximum  likelihood  estimates  (MLE)  of  all  the  unknown  parameters  under  both 
the  hypotheses.  Computing  the  MLE  of  the  signal  amplitude  is  particularly  difficult  in 
the  presence  of  other  unknown  parameters. 

This  paper  proposes  a  Rao  efficient  score  test  [Rao  1948]  as  an  alternative  to  the 
GLRT  for  the  detection  problem  discussed  above.  The  Rao  test  is  shown  to  be  equiva¬ 
lent  to  the  GLRT  when  the  signal  amplitude  is  small.  Consequently,  it  is  also  equivalent 
in  performance  to  a  clairvoyant  Rao  detector,  which  assumes  the  noise  parameters  to  be 
known.  It  greatly  reduces  the  computational  complexity  of  the  GLRT  without  sacrific¬ 
ing  its  optimality  properties.  A  connection  is  established  between  the  Rao  detector  and 
a  locally  optimal  (LO)  detector  which  assumes  the  polarity  of  the  signal  and  the  noise 
parameters  to  be  known.  Computer  simulations  support  the  theoretical  predictions  of 
the  asymptotic  performance  of  the  Rao  detector. 
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The  paper  is  organized  as  follows.  Section  II  introduces  the  Rao  test  as  an  asymp¬ 
totic  approximation  to  the  GLRT.  Section  III  presents  two  different  noise  models  and 
defines  the  GLRT  for  each  one.  Section  IV  derives  the  corresponding  Rao  tests.  Sec¬ 
tion  V  discusses  the  performance  of  the  Rao  test  and  compares  it  to  the  LO  detector. 
Section  VI  derives  the  Rao  detectors  for  the  special  cases  of  autoregressive  (AR)  noise 
with  Gaussian  and  mixed-Gaussian  distributions.  Section  VII  reports  the  results  of 
computer  simulations  and  Section  VIII  sumarizes  the  main  results. 


II.  GLRT  and  Rao  Test 

Let  {ui,u2,---,Ujv}  be  a  set  of  independent  and  identically  distributed  (i.i.d.) 
random  variables  each  having  a  probability  density  function  (PDF)  /(u;  0)  depending 
on  the  q-dimensional  parameter  vector  0  =  [$i  02  •  •  •  0q]T .  Consider  the  following 
composite  hypothesis  testing  problem 


Ho  :  ©  =  [0r  ©arf 

Hi  :  0  =  [0?  0j]r,  0r  #  0 


where  0r  and  0a  are  r  and  s-dimensional  parameter  vectors,  respectively,  with  r+s  =  q. 
0  is  an  r-dimensional  vector  of  0’s.  0a  is  assumed  to  be  unknown.  It  is  sometimes  called 


the  vector  of  nuisance  parameters.  The  GLRT  [Kendall  and  Stuart  1979]  for  testing 
(l)  is  to  decide  Mi  if 


N 

max  TT/K;  ©r,  ©a) 

£<=  = - TT - 


>  1 


max  J]7(un;O,0a) 


n~l 


for  some  threshold  7.  Let  0r  and  0a  be  the  MLE’s  of  0r  and  0a  under  Mi  and  0a  be 
the  MLE  of  0a  under  Mq.  Denoting  the  likelihood  function  by 


N 


T(0r,©a)  =  JJ/(Un;©r,©a), 

71—  1 
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the  statistic  (or  likelihood  ratio )  for  the  GLRT  can  be  written  as 


tG=  £(0rA) 
£(oA) 

The  quantity 

v,(er,e.)  =  Aln£(0rj0<) 

is  defined  as  the  efficient  score  for  the  parameter  and 


(3) 


v(er,e,)  =  (v1(©r,e,)  ^(e,,©,) ...  Fr(©r,0,)]T 

is  the  vector  of  efficient  scores  for  the  parameter  vector  Or*  It  follows  from  the  asymp¬ 
totic  properties  of  the  MLE  that  under  regularity  conditions  of  the  PDF  [Kendall  and 
Stuart  1979]  V(0r,0a)  is  asymptotically  of  the  form 


V(0r,0a)=I(0r,0a)(0-0) 


(4) 


where  ©  is  the  MLE  of  0  and 


[I(0r,0a)]<y=-2? 


d 2 


dBidOj 


In  T(0r,  0a) 


LJ  =  1,2, 


(5) 


which  is  the  ( i,j )  element  of  the  Fisher  information  matrix.  Integration  of  (4)  with 
respect  to  0  results  in  [Cox  and  Hinkley  1974] 


T(0r,0a)  =  Cexp 


—  (0-0)TI(0r,0a)(0-0) 


(6) 


where  C  is  a  constant  not  dependent  on  0.  The  numerator  and  denominator  of  (3) 
correspond  to  the  maximum  of  the  right  hand  side  of  (6)  under  and  Mo,  respectively. 
The  maximization  with  a  few  asymptotic  arguments  yield  [Kendall  and  Stuaxt  1979] 
yeild 

1 


Lg  = 


exp 


-^I©rer(O,0a)0r 


(7) 
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I©r@r  (0r,  0a)  is  the  r  x  r  partition  of  I(0r,  ©a)  corresponding  to  the  parameter  vector 
0r  or 

d2 


I©r©r  (Of}  ©s)  jj'y  =  -E 


<JU  ill  j 


ln£(©r,0s) 


-  1,2,' 


Therefore  for  large  data  records  (asymptotically)  it  follows  from  (7)  that 

2ln£G  =  0^I©r0r(O,0a)0r  (8) 

This  result  was  proved  rigorously  by  Wald  [1943].  An  alternative  approximation  of  iG 
can  be  made  by  rewriting  (6)  the  following  way  by  using  (4) 


T(0r,  0a)  =  Cexp 
The  approximation  of  Iq  is 

lG  =  exp 


ivr(©r,0a)I-1(0r,08)V(0r,0a) 


iv;r(O,0s)J(O,03)Vr(O,0a) 


(9) 


(10) 


where  Vr(0r,  0a)  is  the  r  x  1  partition  of  V(0r,0a)  corresponding  to  0r, 

Vr(0r,0s)  =  [7i(0r,0s)  V2(0r,0s)  ...  Vr(0r,0a)]T 
and  J(0r,  ©a)  is  the  r  x  r  partition  of  I-1(0r,  0a) 

J(©r,©„)  =  [lere,(0r,0.)  -I©,e.(©r,e,)Ii.1©.(©r,6.)lgre.(©r,0.)]"1  (11) 

The  terms  in  the  brackets  of  are  found  by  partitioning  the  Fisher  information  matrix 


for  © 


it©)  = 


I©r©r(Or,©s)  I©r©.  (©r,  ©s) 

I©,  ©r  (Or,  0a)  !©<©>(0r,0a) 


(12) 


and  the  partitions  are  defined  as 

I©r©r(Or,  ©s)  =  E 


/  d  hi  /  \  /  d  In  /  \  T 

\dQr  /  V  dQr  ) 

/31n/\  /ain/\T 

V  <9©r  J  V  dQs  J 


I©r0,(O r,  ©s)  =  E 
I©,©r(Or,©s)=lgr0<(O„Os) 
!©«©,  (Or,  ©s)  —  E 


/ainA/ain/v 

V  5©a  J  V  <9©s  J 


r  x  r 


r  x  s 


s  x  r 


s  x  s 


5 


This  definition  of  the  information  matrices  is  equivalent  to  the  previous  definition  (5) . 
From  (10)  it  follows  that  the  following  test  statistic  proposed  by  Rao  [1948]  is  asymp¬ 
totically  equivalent  to  2  In  to 

tR  =  Vrr(O,03)J(O,0s)Vr(O,0a)  (13) 

Rao’s  test  for  the  complex  hypothesis  testing  problem  (1)  is  to  decide  Mx  if 

Zr  >  i  (14) 

where  qr  is  a  suitable  threshold  and  is  equal  to  2  In  q  if  the  above  asymptotic  equiv¬ 
alence  holds.  The  test  is  called  Rao  efficient  score  test  since  it  uses  the  vector  Vr  of 
efficient  score  functions.  Note  that  estimation  of  0r  and  03  under  Mx  is  avoided  by 
using  the  Rao  test  instead  of  the  GLRT.  This  is  an  outcome  of  the  approximation  of  a 
finite  difference  by  a  derivative  as  in  (4)  which  holds  only  if  the  alternative  hypothesis 
[Hi)  tests  for  small  departures  from  0r  =  0.  The  computational  simplicity  of  the  Rao 
test  makes  it  quite  attractive  for  the  composite  hypothesis  testing  problem  described 
above.  A  more  rigorous  derivation  is  available  in  [Rao  1973]. 

The  statistics  of  Er  are  difficult  to  obtain  in  general.  For  large  data  records  (asymp¬ 
totically)  it  may  be  shown  that  both  2  In  Eq  and  Er  are  distributed  in  the  following 
manner  [Rao  1973]. 

£r  ~  Xr  under  l/0  (15a) 

^■R  ~  X,2(r>  A)  under  Mx  (156) 

Here  x2  represents  a  chi-square  distribution  with  r  degrees  of  freedom  and  x/2(r,  A) 
represents  a  noncentral  chi-square  distribution  with  r  degrees  of  freedom  and  noncen¬ 
trality  parameter  A.  Note  that  x'2(r,0)  =  xl  or  the  distribution  under  M0  is  a  special 
case  of  the  distribution  under  M\  and  occurs  when  A  =  0.  The  noncentrality  parameter 
A,  which  is  a  measure  of  the  discrimination  between  two  hypotheses,  is  given  by 

^  —  ®r  [l©r©r  (0)  ©a)  —  Iere,  (0,  ©a)!©^  (0, 0s)!©,.©,  (0, 08)j  0r  (16) 
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All  the  terms  in  the  brackets  of  (16)  are  partitions  of  I(0r,  0a)  as  given  by  (12). 


III.  The  Detection  Problem  and  the  GLRT  Solution 

III  A.  The  General  Linear  Model 
Consider  the  following  detection  problem. 

l/0:y  =  Wu 

(17) 

l/i:y  =  Wu  +  Ms 

where  s  =  [sx  s2  •••  sN]T  is  a  vector  of  known  signal  amplitudes,  u  =  [ui  u2  •  •  •  uN]T 
is  a  vector  of  i.t.d.  noise  with  a  symmetric  PDF,  /z  is  an  unknown  scalar  (either  positive 
or  negative)  and  W  is  an  invertible  ( N  x  N)  matrix  whose  elements  are  functions  of  a 
set  of  unknown  parameters  ^  tp2  ... 

[W]iy=wy(¥) 

Since  un,  n  —  1, 2,  •  •  •  N  axe  t.t.cf.,  the  PDF  of  u  can  be  expressed  as 

N 

f(u;$)  =  n  /(«*»;*)  (is) 

n=l 

where  f[un  \  $)  is  the  marginal  PDF  of  each  un  dependent  on  the  unknown  parameter 
vector  $.  /  is  assumed  to  be  an  even  PDF,  i.e.,  /(— u)  =  f(u). 

The  linear  model  of  (17)  is  capable  of  representing  a  large  class  of  correlation 
patterns  of  the  background  noise.  The  assumption  of  a  known  PDF  with  unknown 
parameters  $  adds  flexibility  to  the  model  while  still  maintaining  the  parametric  form. 
A  detector  based  on  this  model  would  be  insensitive  to  a  change  of  polarity  of  the  signal 
since  /z  can  be  positive  or  negative.  (17)  is  written  as  the  hypothesis  testing  problem 

U0:OT  =  [0T  ©J] 

:  0T  =  [0T  ©?] 
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Or  7^0 


(19a) 

(196) 


(20) 


where 

Qr  =  n  (a  scalar) 

©s  =  ['F1'  (vector  of  nuisance  parameters) 

The  vector  y  is  a  linear  function  of  the  vector  of  i.i.d.  random  variables  u  under  either 

hypothesis.  Using  this  fact  the  joint  PDF  of  y  is  found  to  be 

N  /  \ 

under  Mo  (21a) 


U»=X^'=  Ww 

n  j  under  M\  (216) 

Un=y~)J-i  u>nj('y)(yj~Pai)  J 


where  wny(^)  are  elments  of  W-1,  which  are  known  functions  of  \F. 

S'W  =  [W-1]ny 

The  GLRT  for  testing  (17)  has  been  shown  to  be  [Kay  and  Sengupta  1986,  3]  equivalent 
to  deciding  Mi  if 

n/f  12u*j  W  te  -  iW ) ;  4 

n=l  \ j  =1 


ta  = 


N  /  N  \ 

I]/ 

n=l  \j= 1  J 


^>1 


(22) 


A  and  $  are  joint  MLE’s  of  ’F,  \i  and  $,  respectively,  under  M\.  $  and  $  are  joint 
MLE’s  of  'F  and  $  under  Mq. 


IIIB.  The  AR  Noise  Model 
The  detection  problem  for  AR  noise  is 

:  y  =  x 

Mi  :  y  =  x  +  /is 


(23) 


with 


x  =  [xi  x2  •  •  •  xN]' 
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It  is  assumed  that  the  sequence  {xi,x2,  •  •  • ,  xN}  is  the  output  of  a  p-th  order  all-pole 
filter  excited  by  white  driving  noise  or 

p 

•^n  ^  ]  ajxn—j  4"  un,  n  =  1, 2,  •  •  • ,  N 

J=1 


Alternately, 


p 

un  ]ajxn-j  j  n  =  1,2,  •  •  •  ,  IV 

j=0 


assuming  ao  =  1.  If  the  samples  y_p+1, y_p+2,  •  •  • , yo  are  assumed  to  be  0,  then  (23) 
can  be  shown  to  be  aspecial  case  of  (7)  where  W”1  is  a  lower  triangular  Toeplitz  matrix 
Kay  and  Sengupta  1986,  3j.  The  results  of  the  previous  section  can  then  be  used  to 
determine  the  GLRT.  Alternatively,  the  GLRT  can  be  derived  on  the  basis  of  the 


conditional  likelihood  function,  assuming  the  data  records  to  be  large.  The  conditional 
likelihood  of  yp+1 , yp+2 ,---,yN  given  yi,y2,---,yp  is 
f  (yp+i  >  yP+ 2 ,  •  •  • ,  yjv  |yi ,  y2,  •  •  • ,  yP) 


N 


II  f  (X/vy»-/;® 

n=p+l  \j=0  , 

n  /(  ~ 

n=p+l  \y=o  J 


under  kb 
under  "Hi 


(24a) 

(246) 


The  GLRT  is  then  given  by 


f-G  = 


N  /  p 

II  f  Yjtj (y«-y  -  hn~j)\ I 

n=p+l  \j=  0 

N  T~v  7 

II  / 


(25) 


n=p+ 1  \y=o  / 

where  hat’s  indicate  MLE’s  under  k0  and  double  hat’s  indicate  MLE’s  under  Hu  re¬ 
spectively.  do  and  do  are  defined  to  be  unity. 

Computation  of  (22)  or  (25)  involves  evaluation  of  the  MLE’s  of  all  the  parameters 
under  H\.  The  Rao  test,  as  indicated  in  the  previous  section,  is  able  to  avoid  this 
computation.  The  following  section  illustrates  how  this  is  accomplished. 
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IV.  Rao  Test  for  the  Detection  Problem 

IVA.  The  General  Linear  Model 

For  the  detection  problem  of  (17)  0r  and  0a  are  given  by  (20).  In  this  case  lR 
defined  in  (13)  can  be  computed  from  the  likelihood  function  or  PDF  of  y  (see  (21)). 
It  is  observed  that 

d lnf(y;  0r,  0a) 
dQr 


Vr(O,0a)  = 


©r=o 


=  ^lnf(y;M>0a) 


/i— 0 


=  |;ln 


d^ln 


N  f  N 

m  d- 

n=l  y=i 
N 


fX= 0 


n=  1 


/i=0 


since  from  (17)  un  can  be  written  as 


N 


u 


n  =  ~  IXSj) 

3-1 


under  Mi.  Therefore 


N 


n- 1  i 


Vr(O,0,)  =  ^ 

n— 1  l 

N  r 

=  £ 

n- 1 


/i— 0 


N  N 

=  E  -!>*•(*)«> 

n=l  \  y=l 


rKj#) 

./K;$) 


(26) 


(27) 


It  has  been  proved  for  the  detection  problem  considered  here  (and  specifically  if  /(u) 
is  an  even  function)  that  [Sengupta  1986]  under  certain  regularity  conditions 


Ier©»  (©r5  ©5)  —  I/x©,  (^j  ©a)  —  0 
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Hence  from  (11)  it  follows  that 


i.e., 


N 

=  T,e 


n=l 

N 


a  (If 


dun 

N 


=  E£ 

n=l 

J-1(O,0,)=/W(O,4',$) 

N  (  N 


-£-■')  (S): 


n=l  y  3=1 


*/(*) 


(28) 


where 


I/($)  =  £ 


ain/\2 

dlin  / 


M=0 


(29) 


Substituting  (27)  and  (28)  in  (13) 


Lr 


N 


E  - 


n=l 


N 

3- 1 


“ 

(  N  *  V 

r 

II 

( s  ^  \ 

/i 

L 

N  N 

E  |  J>(*) 

n=l  \  ;=1 


(30) 


If  $  is  a  scalar,  //(!>)  can  be  stored  in  a  table  as  a  function  of  the  parameter.  But  if 
$  is  a  vector  consisting  of  multiple  parameters,  it  is  a  difficult  task  to  evaluate  it  by 
integration  for  each  value  of  On  the  other  hand  it  is  possible  to  use  the  asymptotic 
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equivalence  [Kay  1985] 


J/($)  =  E 


ft\  2 


N 


N 

E 

n— 1 


( N  .  V 

/'I 

y 

(  N  ^  \ 

/ 

^w»y(^)yy;  $ 

vi=i  y  j 

n  2 


Substituting  in  (30) 


£r  — 


N  (  N 
n=l  \  j  =  i 


y=i 

N 

f  ( 

J- 1 


-1 1  2 


£ 

iV 


N  f  N 
"= 1  \  3= 1 


N 

E 

n=l 


/ 1  £>»,(*)«;  i 

y=1 


-i  2 


In  the  case  of  white  noise 


,  =  /l>  n  =  j,  n,j  =  1,2, ,N 
°>  n,y  =  l,2,...,JV 


Wnj  = 


and  (31)  reduces  to 


Lr  = 


"  r/r(^;$)nn2 


J2Sn 

n= 1 


/(!<»;  *) 


» ( £#)  t 


vn=l  /  n=l 


/'(yn;$) 

/(y»;*) 


as  obtained  by  Kay  [1985]. 


IVB.  The  AR  Noise  Model 


(31) 
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A* 


The  parameter  vector  in  the  AR  noise  case  is 
0r  =  M  (a  scalar) 

©s  =  [aT  <&T]T  (vector  of  nuisance  parameters) 

In  this  case  it  is  easier  to  use  the  conditional  likelihood  functions  given  by  (24).  The 
vector  of  efficient  score  functions  is  found  to  be 


Vr(°,0,)  =  -1„ 


=1'" 


N  p 

II  f  Eay(y»-y-A«»-y);® 

n— p+1  ^y=o 
N 

n  /(«»:*) 

n~p+ 1 


fi= 0 


/i— 0 


where  can  be  written  from  (23)  as 


U*  =  X>(yjW  -  fASn—j) 
3=0 


under  Mi.  Hence 


N 


v,(o.e.)=  £ 

n=p+l 

N 

=  E 

n-p-fl  L 
N 


^  In /(«,;*) 

/  dun \  f  d 


dfj,  J  \dun 

p 


= 

n=p+l  \  yr=0 


p=0 

In /(tin ;  $) 

/  d 


dun 


fi= 0 

ln/(un;$) 


/x=0 


AT 


n=p+l  \  y=0 


fK;$) 

./(«»;*). 


M=0 


(32) 


(33) 


For  the  case  of  AR  noise  it  has  been  proved  that  for  f(u)  an  even  function  and  under 
certain  regularity  conditions  [Sengupta  1986] 


Iere.  (©r,  ©s)  =  Ipe,  (/x,  ©s)  =  0 
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Hence  from  (ll) 


J_1(0r,Oa)  =/MM(M,a,$) 

N 

=  E  E 


n—p+l 

N 

=  E  E 

n=p-fl 


^  dun 
p 


d  ,  ,,  ,A  /  dunx  ' 2 

ln/(un;$)  '  1 


dfj. 


-g-H  (S) 


V 

**e*5  J  (0,  ©5)  =  /MM(0,a,$) 


!/($) 


where  I/($)  is  as  defined  in  (29).  Substituting  (33)  and  (34)  in  (13) 


N 

/  p  \ 2 

E 

~XlQy3rt-j 

n=p-f  1 

V  y 

2 


N 


n=p-fl  \  j—0 


f  (  ^jVn—j]  $ 

J=0 


= 


N  p 

n=p+l  y  j=0  J 


If@) 


(34) 


(35) 


As  before,  evaluation  of  the  integration  involved  in  I j  (4?)  can  be  avoided  by  the  asymp¬ 
totic  equivalence 


I/($)  =  E 


^  N 
iV-p  X^ 

*  n=p+l 


2 


/'  (E8^-* 4 
f  ( Ea>  yn-’]  * 
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Substituting  in  (35) 


£r  = 


N  (  p 

2  ~Y^a>Sn~j 

n=p+l  \  j= 0 


/' 


2 


vz 


1=0 


/  ( Ysm-j-,* 

j=0 


N-p 


N 


'  p  V 

Y?  j  -2>-j 

n=p+l  y  j=0  J 


N 

E 

n=p-f  1 


/' 


2 


^ayy„_j;$ 


VZ 


1=0 


/  |  ^  p-jUn—j'i  ^ 
7=0 


(36) 


Either  (35)  or  (36)  can  be  convenient  for  use  depending  on  the  number  of  unknown 
PDF  parameters. 


V.  Asymptotic  Performance  of  the  Rao  Detector 

Asymptotic  distributions  of  iR  under  M0  and  are  given  by  (15a)  and  (155), 
respectively.  For  the  general  linear  model  0r  =  p  and  08  =  [$r  $T]T  while  for  the 
AR  noise  model  0r  =  p  and  08  =  [a^  Xfoerefore  in  either  case  the  noncentrality 

parameter  is 


^  M  Jnn[0,Qs)  Ip©,  (0,  ©S)I@*@>(0,  0a)I^Q>(O,  0a)] 

=  Al2fpp(0, 0a)  (37) 

since  IM©,  (0, 0a)  =  0,  as  indicated  before.  The  probability  of  false  alarm  is 

Pfa  =  P{Lr  >  V|^o}  (38a) 

and  the  probability  of  detection  is 

Pd  =  P{tR  >  Vl^i}  (385) 
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Both  the  probabilities  can  be  calculated  from  the  tables  of  noncentral  and  central  chi- 
square  distributions,  respectively.  Y  can  be  set  to  produce  a  given  false  alarm  rate  and 
Pd  caii  be  calculated  accordingly. 

It  is  known  that  the  GLRT  for  the  detection  problems  considered  here  is  asymp¬ 
totically  optimal  in  the  sense  that  its  performance  is  equivalent  to  that  of  a  clairvoyant 
GLRT  built  with  perfect  knowledge  of  0S  provided  the  data  record  is  large  [Kay  and 
Sengupta  1986,  3].  This  result  applies  to  the  Rao  test  as  well,  since  it  is  asymptotically 
equivalent  to  the  GLRT.  The  clairvoyant  Rao  test  assumes  that  #  and  $  are  known  or 
that  from  (30) 


lrc  — 


N  N 

y  ( 

n~i  V 


j  N  > 

[ 

n 

\)=1  J 

(N  \ 

r 

2jwny  (*)%•;$ 

i 

^  J 

2 


N  f  N  \ 

n=l  \  j= 1  J 

for  the  general  linear  model  and  in  the  AR  noise  case 


(39) 


,  2 


N  f  p  \ 

Y 

n=p+l  y  j= 0  J 


,;=o 


t-RC  =  — 


\  2 


Nip 

Y  ~YaiSn-j  1  J/($) 

n=p-)-l  \  j~0 


(40) 


Note  that  all  the  parameters  are  assumed  to  be  known.  Asymptotically  tRC  is  dis¬ 
tributed  as 


t-RC  ~  Xr 

under  "Hq 

(41a) 

Lrc  ~  x,2(r,Ac) 

under  "Hi 

(416) 
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where 


Ac  =  0rTI©r©r(O,03)0r 

As  indicated  before,  A  and  Ac  are  asymptotically  equivalent  for  both  the  models  con¬ 
sidered  in  this  paper.  Therefore 


A  Ac  —  fl  7^(0,  ©3) 


(42) 


For  the  general  linear  model  substitution  of  (28)  in  (42)  produces 


A  = 


f  #  A 

2" 

i 

np=l  1 

l  j->  J 

(43) 


while  in  the  AR  case  (35)  can  be  substituted  in  (42)  to  obtain 


A  = 


N 


n=p+l 


y=o 


(44) 


o  is  the  variance  of  un  and  has  been  added  to  facilitate  the  interpretation  of  A.  The 
noncentrality  parameter  is  found  to  dependent  on  the  noise  PDF  only  through  the 
quantity  cr2!/,  assuming  all  the  PDF’s  under  consideration  have  the  same  variance. 
It  is  known  that  this  quantity  is  greater  than  unity  for  all  non-Gaussian  PDF’s  [Kay 
and  Sengupta  1986,  1].  For  a  Gaussian  PDF  <r2 If  =  1.  Therefore  all  other  parameters 


remaining  the  same,  a  non-Gaussian  noise  background  will  produce  a  larger  value  of 
A  than  a  Gaussian  background,  leading  to  a  larger  probability  of  detection.  It  can 
be  shown  that  for  non-Gaussian  and  Gaussian  noise  backgrounds  with  identical  power 
spectral  densities  (PSD),  cr2 If  is  the  ratio  of  the  SNR  necessary  in  the  Gaussian  case 
to  the  SNR  necessary  in  the  Gaussian  case,  in  order  to  produce  a  given  probability 
of  detection  [Kay  and  Sengupta  1986,  3].  Section  VII  reports  the  results  of  computer 
simulations  which  illustrates  this. 


It  is  of  interest  to  notice  the  connection  of  the  Rao  detector  with  a  locally  optimum 
(LO)  detectors  [Middleton  1966],  [Czarnecki  and  Thomas  1984].  A  LO  detector  for  the 
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detection  of  a  known  signal  with  unknown  but  positive  amplitude  in  known  colored  noise 
maximizes  the  probability  of  detection  in  the  neighborhood  of  p  =  0  and  is  represented 
by  the  statistic 

Tlo  =  Yim—  Inf  (y;ju,0a)  (45) 

with  usual  notations.  The  test  decides  Mi  if 

Tlo  >  Ilo 

where  7 lo  is  an  appropriate  threshold.  An  equivalent  test  is  to  decide  Mx  if 

^(0,  ©.) =  t'lo  > = 11  A'‘(0’  ®*)  <46) 

Since  the  performance  of  the  Rao  detector  is  asymptotically  equivalent  to  a  correspond¬ 
ing  clairvoyant  detector,  it  is  possible  to  do  a  comparison  on  the  basis  of  known  nuisance 
parameters.  From  (13),  tRC  is  given  by 

Zrc  =  [VrT  (0, 0a)  J(0,  ©,)' Vr  (0,  ©,)] 

~  d  1 2 

—  lnf(y;0,©8) 

05) 

which  is  observed  to  be  the  square  of  T!lo.  Tlo  is  Gaussian  for  large  data  records 
[Middleton  1966]  while  Zrc  is  chi-square  distributed  with  a  single  degree  of  freedom. 
Therefore  for  a  given  probability  of  false  alarm  the  locally  optimal  detector  will  have 
a  slightly  larger  probability  of  detecting  small  amplitudes  of  signal.  This  makes  intu- 
tive  sense,  because  (46)  is  a  one-sided  test  while  GLRT  and  the  Rao  detectors  allow 
for  positive  and  negative  amplitudes  of  the  signal.  This  is  an  example  of  trading  off 
performance  in  order  to  make  the  test  two-sided.  This  comparison  applies  to  the  case 
of  known  nuisance  parameters  only.  In  the  case  of  unknown  nuisance  parameters  the 
LO  detector  is  not  defined,  while  the  Rao  detector  exists  and  achieves  the  performance 
of  a  clairvoyant  detector  asymptotically.  Furthermore,  the  LO  detectors  for  positive 
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and  negative  values  of  n  are  different  (i.e.,  the  inequality  is  reversed  for  negative  am¬ 
plitudes)  and  hence  they  are  impractical  in  situations  where  the  polarity  of  the  signal 
can  change. 

VI.  Rao  Test  for  Gaussian  and  Mixed-Gaussian  Noise 

VIA.  The  Gaussian  Case 


The  Rao  detector  is  now  derived  for  the  special  case  of  the  zero  mean  Gaussian 
noise.  Attention  is  restricted  to  AR  noise  only  because  of  the  abundance  of  available 
results  [Bowyer  1979],  [Kay  1983].  In  this  case  $  =  cr2,  a  scalar,  and 


/(^n 


such  that 


f(«njg2)  _  Un 


/(«.;* 2)  - 

Substituting  (47)  in  (35)  and  using  the  fact  that  If  {a2)  =  ifo 2  for  a  Gaussian  PDF, 

(  Sn  +  I  I  l In  + 


f-R  = 


n=p+l 


y=i 


N 


(48) 


a  (  Sn  +  y^dySn-y 


n=p+l 


j-l 


where 

1  N 

®  =  jy  _  „  (Vn  +  ajVn-f)2 

Fn=p+ 1 

and  the  estimates  are  obtained  under  Mo  or  assuming  the  signal  to  be  absent. 

Figure  1  shows  the  block  diagram  of  the  Rao  detector.  In  the  Gaussian  case  the 
least  squares  estimators  are  known  to  be  close  to  the  MLE  and  hence  any  of  the  least 
squares  techniques  (such  as  the  autocorrelation  method,  the  covariance  method,  the 
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Forward/backward  method  and  so  on)  [Kay  1986]  can  be  used  to  produce  the  estima¬ 
tors  in  (48).  The  numerator  appears  to  be  a  prewhitener-correlator  [Van  Trees  1968] 
except  for  the  squaring  required  to  make  the  test  two-sided.  The  denominator  serves 
to  normalize  the  statistic  so  that  a  constant  false  alarm  rate  (CFAR)  is  maintained 
asymptotically.  Note  that  it  is  not  necessary  to  estimate  p  or  to  assume  a  prior  value 
for  it.  Also,  assuming  ai,  <12,  •  •  • ,  ap  =  0  or  white  noise  will  result  in  the  square  a  sim¬ 
ple  correlator  or  a  matched  filter  in  the  numerator  of  lR,  which  is  known  to  be  the 
uniformly  most  powerful  (UMP)  test  for  the  detection  of  a  known  signal  in  white  noise 
[Van  Trees  1968].  In  the  case  of  the  general  linear  model  the  moving  average  (MA) 
filter  with  transfer  function  A(Z )  defined  by 

p 

A(Z)  =  l  +  J2a:Z~j 
3= 1 

will  have  to  be  replaced  by  an  inverse  transformation  W-1. 

VIB.  The  Mixed- Gaussian  Case 


The  mixed-Gaussian  PDF  is  particularly  useful  in  representing  a  special  class  of 
non-Gaussian  noise  processes,  namely,  a  nominally  Gaussian  background  contaminated 
with  occasional  impulses  [Sengupta  and  Kay  1986,  lj.  The  PDF  is  given  by 

f  4 

/(«n;g|,g/,g)  =  -^L^eV  2aB 

najj 

where  e  is  a  mixture  parameter  and  0  <  e  <  1.  The  subscripts  B  and  I  are  used  to 
denote  background  and  interference,  respectively.  An  alternative  formulation  of  the 


PDF  is 


ut 


u : 


/(un;«r|,/9,e)  = 


1  —  e 


V2 


2ob)  + 


7TOi 


\/27r  po% 


2Pa% 


(49) 
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where  p  —  crj/aB.  Assuming  p  »  1,  one  can  allow  for  a  wide  range  of  amplitudes  and 
frequencies  of  occurence  of  spikes  by  appropriately  choosing  p  and  e.  The  Rao  detector 
proposed  in  this  paper  is  suitable  for  the  mixed-Gaussian  PDF  (since  it  is  symmetric) 
and  therefore  is  applicable  to  many  sonar  and  radar  detection  problems  where  the  noise 
background  contains  impulses  resulting  from  reverberation  and  clutter,  respectively 
[Kay  1983].  The  Rao  detector  for  mixed-Gaussian  PDF  is  now  derived.  Once  again 
only  the  AR  noise  case  is  considered.  It  is  assumed  that  <j|j  and  p  are  known,  so  that 
$  =  e,  a  scalar.  From  (49)  it  follows  that 


/'K*; e) 

/(««;  e) 


Ur 


ut 


1-e 


V2 


7rcri 


2°b)  + 


p^j2-Kpo2B 


1-e 


K 

2 


+ 


un  tv  \ 

--3-r(un;e) 


V27rPaB 


ut 


2  pa 


B , 


ut 


2Pgb 


(50) 


where  r(un;e)  represents  the  term  in  the  brackets.  Note  that  (50)  differs  from  (47) 
only  by  the  weighting  function  T(un;e)  which  suppresses  the  high  values  of  un  [Kay 
and  Sengupta  1986,  2].  In  the  Gaussian  case  e  =  0  and  T(un;0)  =  1  for  all  un.  The 
weighting  function  T(un;e),  which  is  an  even  function  of  un,  thus  accounts  for  the 
non-Gaussian  nature  of  the  PDF.  In  the  mixed-Gaussian  case,  lR  can  be  calculated  by 
substituting  (50)  into  (35) 


t-R  = 


N 

X^  (  Sn  +  X^  I  [  Vn  +  \  T  (  yn  +  ayJ/n-j ;  6 


p 


1  2 


3= 1 


N 


aB^f(^)  ^  (  Sn  +  sn-j 


n=p+l 


3=1 


(51) 

Note  that  T  and  1/  depend  on  the  value  of  e.  Figure  8.2  shows  the  block  diagram 
corresponding  to  (51)  and  has  two  apparent  differences  with  Figure  1  (the  Gaussian 
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detector).  These  are  the  blocks  used  to  compute  T  and  If.  Substitution  of  cr|,  T  and 
If  by  a2,  1  and  1/d2,  respectively,  would  lead  one  back  to  the  Gaussian  case  described 
by  (48)  and  Figure  1.  The  function  F  can  be  replaced  by  a  suitable  approximation, 
such  as  the  Butterworth function  to  avoid  computation  of  the  exponentials.  The  MLE’s 
can  also  be  replaced  by  a  reasonably  accurate  estimator,  (e.gr.,  a  weighted  least  squares 
estimator)  [Kay  and  Sengupta  1986,  2].  With  all  these  simplifications  (51)  would  pro¬ 
duce  a  detector  which  will  reduce  computation  substantially  without  a  significant  loss 
of  performance. 

VTI.  Computer  Simulations  of  Performance  of  the  Rao  Detector 

This  section  reports  the  results  of  computer  simulations  of  the  performance  of  the 
Rao  detector.  The  clairvoyant  Rao  detector  is  considered  to  be  one  basis  of  comparison, 
while  the  theoretical  or  asymptotic  performance  is  regarded  as  another.  Two  AR  noise 
processes  are  selected  for  computer  simulations.  The  corresponding  parameters  are 
listed  in  Table  A.  Process  I  is  broadband  while  process  II  is  narrowband.  A  mixed- 
Gaussian  process  with  <t|j  =  1,  p  =  100  and  e  =  0.1  is  chosen  as  the  driving  noise,  e  is 
assumed  to  be  unknown  and  so  estimated,  while  <7q  and  p  are  assumed  to  be  known. 
The  known  part  of  the  signal,  t.e.,  sn,  n  =  1,2,  •  •  •  ,N  is  assumed  to  be  unity,  p 
is  adjusted  to  yield  different  values  of  SNR  in  the  following  way.  The  theoretical  noise 
variance  is,  as  obtained  from  (49), 

a2  =  a2B[(l  -  e)  +  ep]  (52) 

The  noise  power  Pn  of  the  AR  process  can  be  obtained  from  the  stepdown  procedure 
using  o2  and  the  process  parameters  [Kay  1986].  The  signal  power  (actually  the  signal 
energy)  is  the  defined  as 

Ps  =  Np2 

Defining  SNR  to  be  the  ratio  of  Ps  and  Pn  (which  represents  the  SNR  at  the  output 
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of  a  correlator) ,  it  follows  that 


/  (SNR)Pn 

»  =  n  (53) 

Thus  fx  is  calculated  for  a  given  process  such  as  to  produce  a  desired  SNR.  N,  the 
number  of  data  points  is  chosen  to  be  1000.  A  probability  of  false  alarm  PFa  =  0.01  is 
used  to  evaluate  the  detection  performance.  The  value  of  7'  necessary  for  this  is  6.635, 
as  obtained  by  a  search  routine  so  as  to  satisfy  PFa  =  0.01. 

In  order  to  evaluate  the  asymptotic  performance  of  the  Rao  test,  the  noncentrality 
parameter  A  =  Ac  as  given  by  (44)  is  calculated  from  fx  for  each  value  of  SNR  as  per 
(53).  For  the  chosen  values  of  the  PDF  parameters,  o2If  is  calculated  by  numerical 
integration  and  is  found  to  be  9.0.  PD  as  defined  by  (386)  is  computed  from  a  table 
using  the  values  of  A  and  7'  obtained  previously. 

The  theoretical  value  of  the  threshold  7',  as  described  above,  is  used  for  the  Rao 
detector.  This  requires  one  to  verify  of  the  theoretical  predictions  of  the  asymptotic 
statistics  of  under  )lo  in  order  to  use  the  theoretical  threshold.  Computer  simulation 
results  based  on  1000  experiments  with  1000  data  points  each  using  the  theoretical 
threshold  in  the  absence  of  the  signal  result  in  a  false  alarm  in  11  cases  for  process  I 
and  10  cases  for  process  II,  corresponding  to  an  experimental  false  alarm  rate  of  0.011 
and  0.01,  respectively,  which  are  very  close  to  the  true  value  of  0.01.  This  is  expected, 
since  the  MLE’s  when  H0  is  true  are  expected  to  be  more  accurate  than  the  MLE’s 
when  is  true  (since  the  MLE’s  are  computed  under  the  assumption  that  V,Q  is  true). 
Therefore  the  use  of  a  theoretical  threshold  is  justified  and  convenient.  The  statistic 
£r  is  computed  from  (51)  for  500  different  blocks  of  data,  each  of  length  N  =  1000, 
for  a  given  SNR.  The  number  of  times  the  statistic  exceeds  7',  scaled  by  500  (the 
number  of  experiments),  is  regarded  as  the  experimental  value  of  the  probability  of 
detection.  This  is  repeated  for  different  values  of  SNR  in  a  suitable  range  (so  as  to 
observe  the  transition  from  PD  =  0  to  PD  =  1).  The  MLE  of  the  AR  filter  parameters 


23 


involved  in  (51)  under  "Hq  are  replaced  by  the  two-stage  weighted  least  squares  estimator 
proposed  by  the  authors  [Kay  and  Sengupta  1986,  2]  to  reduce  computation  and  avoid 
convergence  problems  for  short  data  records.  Once  these  parameters  are  estimated ,  e 
is  computed  as 


from  the  prediction  error  power  d2  of  an  MA  filter  fed  by  the  observed  data  with 
coefficients  &i,  a2,  •  •  •  ap  (see  (52)).  I/(e)  is  computed  for  each  e  from  a  stored  table  by 
interpolation.  T  in  the  numerator  of  tR  is  replaced  by  the  Butterworth  approximation 


f(u)  — - - — g  +  k2 

u 


1+  — 
uc 


in  order  to  reduce  computation.  The  parameters  Ku  K2,  uc  and  (3  are  selected  to 
approximate  T  over  the  range  of  e. 

The  statistic  lRC  for  the  clairvoyant  Rao  test,  as  defined  by  (40),  is  computed 
in  a  similar  way  using  true  values  of  all  the  parameters.  This  is  used  to  determine 
the  probability  of  detection  ( PD )  of  the  clairvoyant  Rao  detector  for  each  SNR  via 
computer  simulations. 

Figure  3  plots  the  probabilities  of  detection  of  the  Rao  detector  and  the  clairvoyant 
Rao  detector  along  with  the  theoretical  or  asymptotic  performance  vs.  SNR  for  the 
noise  process  I.  Figure  4  plots  the  same  for  the  noise  process  II.  The  three  curves 
exhibit  remarkable  similarity  for  each  process.  It  can  be  concluded  that  (15)  and 
(38)  adequately  represent  the  performance  of  the  computationally  efficient  detector 
proposed  in  this  paper,  all  the  approximations  incorporated,  even  for  moderately  sized 
data  records. 

The  performance  of  the  Rao  detector  is  now  compared  to  that  of  a  Gaussian  Rao 
detector.  A  Gaussian  Rao  detector  (defined  by  (48))  is  computationally  much  simpler. 
The  purpose  of  comparing  these  two  detectors  is  to  show  how  much  one  can  lose  as  a 
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result  of  an  incorrect  Gaussian  assumption  about  the  driving  noise.  All  least  square 
estimators  are  close  approximations  to  the  MLE  if  the  PDF  is  assumed  to  be  Gaussian. 
The  covariance  method  is  used  to  obtain  the  estimates  involved  in  (48) .  The  same 
theoretical  threshold  7'  is  used.  (Computer  simulations  based  on  1000  experiments 
with  1000  data  points  each  using  the  theoretical  threshold  indicate  for  process  I  an 
experimental  false  alarm  rate  of  0.009  which  is  close  to  the  true  value  of  PFA ).  The 
experimental  threshold  for  process  II  is  found  to  be  0.012  PD  is  calculated  from  the 
number  of  times  (out  of  500  independent  trials)  iR  exceeds  7'.  N  =  1000  is  used 
and  PD  is  computed  for  different  values  of  SNR  in  a  suitable  range.  Figures  5  and 
6  plot  the  resultant  experimental  performance  of  the  Gaussian  Rao  detector  in  the 
same  scale  with  those  of  the  mixed-Gaussian  Rao  detectors,  actual  and  asymptotic, 
for  the  two  noise  processes  described  before.  They  indicate  a  substantial  degradation 
in  performance  as  a  result  of  the  Gaussian  assumption  of  the  driving  noise.  It  is  also 
observed  that  the  performance  of  the  Gaussian  Rao  detector  matches  the  predicted 
asymptotic  performance  of  a  Gaussian  Rao  detector  in  Gaussian  noise  having  equivalent 
variance.  This  can  be  explained  in  the  following  way.  Since  it  is  assumed  that  sn  = 
1,  n  =  1, 2,  •  •  • ,  N  (48)  can  be  written  as 


(54) 


Assuming  the  estimates  of  the  AR  filter  parameters  to  be  reasonably  accurate,  the 
term  in  the  parantheses  is  approximately  the  nth  sample  of  the  driving  noise  which  is 
assumed  in  this  case  to  be  mixed-Gaussian.  By  Central  Limit  Theorem  arguments  the 
sum  of  N  -  p  such  random  variables  scaled  by  (N  -  p)d2  has  a  Gaussian  PDF  with 
variance  one  assuming  d2  to  be  close  to  the  variance  of  the  said  mixed-Gaussian  PDF. 
Therefore  the  asymptotic  statistics  of  tR  in  this  case  will  be  the  same  as  what  it  would 
have  been  if  the  noise  were  Gaussian  with  the  same  variance.  It  can  be  concluded 
that  the  asymptotic  statistics  of  the  Gaussian  Rao  detector  in  mixed-Gaussian  noise  is 


25 


represented  by  (15)  and  (44),  despite  modeling  error.  However  this  conclusion  holds 
for  large  data  records  and  small  D.C.  level  signal  amplitudes  only.  A  large  value  of 
is  expected  to  cause  a  significant  degradation  of  performance  of  the  LS  estimators  so 
that  the  above  interpretation  of  (54)  is  no  longer  valid.  A  short  data  record  will  make 
the  Central  Limit  Theorem  inapplicable.  Comparing  the  asymptotic  performances  of 
the  Gaussian  and  non-Gaussian  detectors,  o 2  If  is  found  to  describe  quantitavely  the 
improvement  of  the  mixed-Gaussian  detector  over  the  Gaussian  detector  when  the  true 
noise  actually  fits  the  mixed-Gaussian  model  and  the  signal  is  a  D.C.  level.  Figures  5 
and  6  show  a  constant  difference  of  approximately  10  dB  between  the  performances  of 
the  Gaussian  and  mixed-Gaussian  detectors  which  matches  the  theoretical  prediction 
of  9.6  dB  (since  a2 If  «  9  in  this  case). 

VIII.  Summary 

The  Rao  efficient  score  test  proposed  in  this  paper  is  found  to  be  well  suited  for 
the  problem  of  detecting  a  weak  signal  of  unknown  amplitude  in  the  presence  of  col¬ 
ored  non-Gaussian  noise  of  unknown  PDF  and  PSD  parameters.  Since  the  Rao  test  is 
asymptotically  equivalent  to  the  GLRT,  it  shares  all  the  attractive  asymptotic  prop¬ 
erties  possessed  by  the  GLRT.  It  greatly  reduces  computation  by  completely  avoiding 
estimation  of  the  unknown  parameters  under  The  Rao  detectors  were  derived  for 
Gaussian  and  mixed-Gaussian  background  noise  processes.  The  performance  of  the 
Rao  detector  is  found  to  be  equivalent  to  that  of  a  clairvoyant  Rao  detector  built  with 
perfect  knowledge  of  the  nuisance  parameters.  The  experimental  performance  matches 
the  theoretical  predictions  of  the  asymptotic  performance  for  large  data  records.  A 
detector  which  assumes  the  noise  PDF  to  be  Gaussian,  while  it  is  actually  not  so,  is 
found  to  be  much  inferior  in  performance. 
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Table  A:  Parameters  of  the  AR  processes  used  for  simulation 


Process 

ai 

C&2 

<*3 

a4 

poles 

i 

-1.352 

1.338 

-0.662 

0.240 

0.7  exp[y27r(0.12)] 

0.7  exp[y27r(0.2l)] 

ii 

-2.760 

3.809 

-2.654 

0.924 

0.98exp[y27r(0.11)j 

0.98exp[y2?r(0.14)] 

noise 


prewhitener 


correlator 


numerator 


Figure  1 


Block  diagram  of  Rao  detector  for  Gaussian  noise 
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Figure  3  Performance  of  the  Rao  detector  for  noise  process  I 
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Figure  4  Performance  of  the  Rao  detector  for  noise  process  II 
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Figure  5 


Performance  of  the  Gaussian  and  mixed-Gaussian 
Rao  detectors  for  mixed  Gaussian  noise  process  I 
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Figure  6 


Performance  of  the  Gaussian  and  mixed-Gaussian 
Rao  detectors  for  mixed-Gaussian  noise  process  II 
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Figure  6  Performance  of  the  Gaussian  and  mixed-Gaussian 

Rao  detectors  for  mixed-Gaussian  noise  process  II 
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Figure  3  Performance  of  the  Rao  detector  for  noise  process  I 


OFFICE  OF  NAVAL  RESEARCH 
STATISTICS  AND  PROBABILITY  PROGRAM 

BASIC  DISTRIBUTION  LIST 
FOR 

UNCLASSIFIED  TECHNICAL  REPORTS 


FEBRUARY  1982 


Copies 


Statistics  and  Probability 
Program  (Code  41 1 (SP ) ) 

Office  of  Naval  Research 
Arlington,  VA  22217  3 

Defense  Technical  Information 
Center 

Cameron  Station 

Alexandria,  VA  22314  12 

Commanding  Officer 
Office  of  Naval  Research 
Eastern/Central  Regional  Office 
Attn:  Director  for  Science 
Barnes  Building 
495  Summer  Street 
Boston,  MA  02210  1 

Commanding  Officer 
Office  of  Naval  Research 
Western  Regional  Office 
Attn:  Dr.  Richard  Lau 
1030  East  Green  Street 
Pasadena,  CA  91101  1 

U.  S.  ONR  Liaison  Office  -  Far  East 
Attn:  Scientific  Director 
APO  San  Francisco  96503  1 

Applied  Mathematics  Laboratory 
David  Taylor  Naval  Ship  Research 
and  Development  Center 
Attn:  Mr.  G.  H.  Gleissner 
Bethesda,  Maryland  20084  1 

Comnandant  of  the  Marine  Coprs 
(Code  AX) 

Attn:  Dr.  A.  L.  Slafkosky 
Scientific  Advisor 
Washington,  DC  20380  1 


Copies 


Navy  Library 

National  Space  Technology  Laboratory 

Attn:  Navy  Librarian 

Bay  St.  Louis,  MS  39522  1 

U.  S.  Army  Research  Office 
P.0.  Box  12211 
Attn:  Dr.  J.  Chandra 
Research  Triangle  Park,  NC 

27706  1 

Di rector 

National  Security  Agency 

Attn:  R51 ,  Dr.  Maar 

Fort  Meade,  MD  20755  1 

ATAA-SL,  Library 
U.S.  Army  TRADOC  Systems 
Analysis  Activity 
Department  of  the  Army 
White  Sands  Missile  Range,  NM 
88002  1 

ARI  Field  Unit-USAREUR 

Attn:  Library 

c/o  ODCSPER 

HQ  USAEREUR  &  7th  Army 

APO  New  York  09403  1 

Library,  Code  1424 
Naval  Postgraduate  School 
Monterey,  CA  93940  1 

Technical  Information  Division 
Naval  Research  Laboratory 
Washington,  DC  20375  1 

OASD  (I&L),  Pentagon 

Attn:  Mr.  Charles  S.  Smith 

Washington,  DC  20301  1 


Copies  Copies 


Di rector 
AMSAA 

Attn:  DRXSY-MP,  H.  Cohen 
Aberdeen  Proving  Ground,  MD  1 
21005 

Dr.  Gerhard  Heiche 
Naval  Air  Systems  Comnand 
(NAIR  03) 

Jefferson  Plaza  No.  1 
Arlington,  VA  20360  1 

Dr.  Barbara  Bailar 
Associate  Director,  Statistical 
Standards 
Bureau  of  Census 
Washington,  DC  20233  1 

Leon  Slavin 

Naval  Sea  Systems  Command 
(NSEA  05H) 

Crystal  Mali  #4,  Rm.  129 
Washington,  DC  20036  1 

B.  E.  Clark 

RR  #2,  Box  647-B 

Graham,  NC  27253  1 

Naval  Underwater  Systems  Center 
Attn:  Dr.  Derrill  J.  Bordelon 
Code  601 

Newport,  Rhode  Island  02840  1 

Naval  Coastal  Systems  Center 
Code  741 

Attn:  Mr.  C.  M.  Bennett 
Panama  City,  FL  32401  1 

Naval  Electronic  Systems  Command 
(NELEX  612) 

Attn:  John  Schuster 
National  Center  No.  1 
Arlington,  VA  20360  1 

Defense  Logistics  Studies 
Information  Exchange 
Army  Logistics  Management  Center 
Attn:  Mr.  J.  Dowling 
Fort  Lee,  VA  23801  1 


Reliability  Analysis  Center  (RAC) 

RADC/RBRAC 

Attn:  I.  L.  Krulac 

Data  Coordinator/ 

Government  Programs 
Griff iss  AFB,  New  York  13441  1 

Technical  Library 

Naval  Ordnance  Station 

Indian  Head,  MD  20640  1 

Library 

Naval  Ocean  Systems  Center 

San  Diego,  CA  92152  1 

Technical  Library 
Bureau  of  Naval  Personnel 
Department  of  the  Navy 
Washington,  DC  20370  1 

Mr.  Dan  Leonard 
Code  8105 

Naval  Ocean  Systems  Center 

San  Diego,  CA  92152  1 

Dr.  Alan  F.  Petty 
Code  7930 

Naval  Research  Laboratory 
Washington,  DC  20375  1 

Dr.  M.  J.  Fischer 
Defense  Communications  Agency 
Defense  Communications  Engineering 
Center 

1860  Wiehle  Avenue 

Reston,  VA  22090  1 

Mr.  Jim  Gates 
Code  9211 

Fleet  Material  Support  Office 
U.  S.  Navy  Supply  Center 
Mechanicsburg,  PA  17055  1 

Mr.  Ted  Tupper 
Code  M-311C 

Military  Sealift  Command 
Department  of  the  Navy 
Washington,  DC  20390  1 


